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The viscoelastic effect on the hydrodynamic relaxation in semidilute polymer solutions is 
investigated. From the linearized two-fluid model equations, we predict that the dynamical 
asymmetry coupling between the velocity fluctuations and the viscoelastic stress influences on 
the hydrodynamic relaxation process, resulting in a wave-number-dependent shear viscosity. 
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The viscoelasticity in entangled polymer solutions has been extensively studied from various 
viewpoints.BS When we investigate the mesoscopic dynamics of polymer solutions, the two- fluid 
modeli'BB gives a good starting point. The two-fluid model was put forward by Brochard and de 
Gennes.i'& Later, the current form of the two-fluid model was derived by Doi and Onuki.& Its 
validity has been recognized through researches of a number of intriguing viscoelastic phenomena 
such as shear- induced concentration fluctuationsJli'B& nonexponential decay in dynamic scatter- 
ing near equilibrium,!' i'& phase separation dynamics, and so on. The two-fluid model 
contains the following important concepts for understanding the above various phenomena, (i) In 
entangled polymer solutions, the network stress can act on the two components asymmetrically. 
That is, in polymer solutions, the stress is supported by the polymer chains much more than by 
the solvent, and the resultant imbalance of the stress causes the relative motion of the two compo- 
nents. This effect is called dynamical asymmetry coupling.!' Bi^ (ii) Considering the viscoelastic 
effect on the transport phenomena in polymer solutions, the cooperative motion of the network 
structure resulting from the entanglement of the polymers is also important. Such a cooperative 
scale intrinsic to entangled polymer solutions is characterized by the so-called viscoelastic length 
S,ve-, which was first introduced by Brochard and de Gennes.i'B 

Previously, Doi and Onukii showed that the diffusion process is drastically influenced by the 
dynamical asymmetry coupling between the concentration fluctuations and the longitudinal mode 
of the viscoelastic force (stress-diffusion coupling). Their analysis predicted an additional wave 
number dependence of the kinetic coefficient for concentration due to the viscoelasticity, which was 
confirmed by recent experimental studies (see ref. 14 for example). However, in their analysis, very 
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viscous systems without the velocity fluctuations were assumed, so the effect of the viscoelasticity 
on the relaxation of the velocity fluctuations (hydrodynamic relaxation) was not investigated. Here 
we take the transverse parts of the viscoelastic force and the velocity fluctuations into consideration. 
We shall then show that the dynamical asymmetry coupling between these transverse modes has 
a strong influence on the hydrodynamic relaxation, resulting in a wave-number-dependent shear 
viscosity. 



In the present analysis we use the two-fluid modeld'OLF as the basic equations describing the 
dynamics of polymer solutions. Here we briefly survey the two-fluid model equations. The readers 
who want to know the details of the two-fluid model, see the original papei^^ or review articles. El) 

Let Vp(r, t) and Vs(r, t) be the average velocities of polymer and solvent, respectively, and tp{r, t) 
be the volume fraction of polymer at point r and time t. ip satisfies a conservation law with v^: 

^V = -V.(^v,). (1) 

The volume- averaged velocity v = ipVp + (1 — 'ip)vs obeys the following hydrodynamic equation 

(9v dF 

= -Vp + r/oV^v-^V— + V-aP, (2) 

where p is the average mass density, r]Q is the solvent viscosity, F is the free-energy functional, 
(T^ is the viscoelastic stress, and the hydrostatic pressure p is determined by the incompressibility 
condition 

V • V = 0. (3) 

For slow motions, the two-fluid model gives the convective velocity of polymer Vp as 

where C(V') is the friction coefficient between the two components, and is given by 

C(^) = e^r^oC'. (5) 

Here is the blob length of order b/ip, with b being the monomer size. Equations (1) and (4) 
imply that the nonvanishing viscoelastic stress produces a diffusion. This effect is called the stress- 
diffusion coupling. The viscoelastic stress due to the network deformations is given by 

aP = G{ij)W -{W -6). (6) 

This form is derived from the free-energy functional introduced in eq. (8) shown below. Here 
G{ip) is the V'-dependent shear modulus, S denotes the unit matrix and W is interpreted as being 
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a long-lived strain variable. We assume that the equation of motion of W is given by the following 
dynamic equation, EE) 

-Qf + i^p- = (Vvt .w + W■V^^p)- -{W - 6), (7) 
where T^ip) is the V'-dependent relaxation time of shear stress. 

Now we shall investigate the relaxation dynamics of concentration fluctuations and velocity fluc- 
tuations within the linearized approximation of the above set of equations. We set cq = c{^q), 
To = t(^o) and Go = G{ipo) with ipQ being the average concentration. The free-energy functional 
is assumed to be of the following Ginzburg-Landau type 

F{i;,W}= I dr[Ko^^ + ^co\VM^ + \GoY.5W^^], (8) 

ij 

where 6Wij = Wij — 5ij. The Fourier transform of an arbitary function g{r) is defined by 

<7k= / dre-^%(r). (9) 



In terms of the Fourier components, the linearized equations are given by 

= + 1(5 - kk) • k ■ al (10) 

d 

= -rk'^V'k + OiLk • k • <T^, (11) 
d 1 

+2Goarkfei%(5V'k - Goa^Lk^ E(^'^'<,/j + ^^^^ 
= GoSWi,^^^, (13) 
where vj^ is the transverse part of Vk, k = k/Zc, and 

« = ^ (14) 
Wo 

is the dynamical asymmetry parameter of the polymer solutions. The kinetic coefficient L is 
given by 

L = Vg(l-Vo)VC(V'o). (15) 
The decay rate in the absence of the viscoelastic coupling is given by 

rk = L(ro + cofc2)/c2. (16) 

The viscoelastic force f is defined as 

f = V • (jP. (17) 
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For the present purpose, it is convenient to decompose the Fourier transform fj^ of the viscoelastic 
force into longitudinal and transverse components, 



fk = kkfk + (5-kk)fk. 



(18) 



Let us set 



fl^=(^--kk)fk. 



(19) 
(20) 



From the linearized eqs. (10)- (13), we obtain the following two sets of equations of motion for the 
longitudinal modes and the transverse modes. 
For the longitudinal modes, we obtain 



d 

— (5V'k = -rk#k - aLZk, 

|-Zk = --(1 + 2e|,fc2)Zk - 2Goar^kHi;^, 
dt tq 

where Zj^ = f|j and ^^e is a viscoelastic length or a magic length defined by 



(21) 
(22) 

(23) 



This length was first introduced by Brochard and de Gennes.Bo' The physical significance of the 
viscoelastic length has been discussed by many researchers (see refs. 10 and 13 for example). 
For the transverse modes, we obtain 



dt 



T]ok 



1 



Vk + -fk > 

p p 



|fk = + ^'ek')f^ - Gok'^ri. 



(24) 
(25) 



From eqs. (21) and (22), the dynamic structure factor S]^{t) is given by 



Skit) = < <5Vk(i)'JV'-k(0) > 



^k(O) 

10+ — LO- 



To 



To 



— UJ—t 



(26) 



where < • • • > denotes the equilibrium average and 
1 



UJ± 



1 



To 



Fk + -(1 + 2^^2)1 '-4^^ 

To J Tq 



(27) 



This form of the dynamic structure factor was first proposed by Brochard and de GennesB& on 
the basis of a phenomenological argument. Afterwards, Doi and OnukiB* derived the same form of 
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the structure factor by analysing their two-fluid model equations which is employed in the present 
analysis (see also ref. 11). 



The viscoelasticity is characterized by the complex shear modulus G^{uj) = G'y,{uj) + 
where G'-^^{lo) and G'-^{io) are called the storage modulus and the loss modulus, respectively .11* Since 
eqs. (24) and (25) are linear, we may assume, without loss of generality, that vj^(t) and (^{t) 
depend on time as 



vi^(i) = vi^e-*, (28) 

From eqs. (24) and (25), the total shear stress o"k(t) is expressed as 

iuj 

where [kvj^^]''' is the transpose of the dyadic kvj^^. In the present system, GJ^(a') and G'-I^{uj) are 
given by 



fk (0 = fk e-*. (29) 



^k(i) = + i[kv^]^}e-*, (30) 



G'^(c.) = u; [^0 + Goto Wt-f^ek^r 1 ' ^''^ 

which are non-local in both space and time. The complex shear viscosity is given by 

( ^ ^k(^) ,oo, 
Vki^) = — • 33 

The imaginary part of G^(lj) gives the shear viscosity for long-time scale motion (u = 0) as 

77k(0) = lim = ^0 + J (34) 

where we have defined rjp = GqTq. The asymptotic behavior of ^(0) is represented by 

VkiO) = Vo + Vp (k^ve^l), (35) 

For polymer solutions, the two-fluid model gives the viscoelastic length a 

^ve = {r]p/r]o)kb- (37) 

Thus, eq. (36) is also expressed as 

r?k(0) = r/o(l + ^), (A:e.e»l), (38) 
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which exphcitly shows that the hydrodynamic interaction becomes weak beyond the blob length 
^fe. It is worth mentioning that this property is due to the dynamical asymmetry coupling (a ^ 0) 
between the viscoelastic stress and the velocity fluctuations. For polymer solutions, it is well known 
that the difl'usion process is drastically influenced by the dynamical asymmetry coupling between 
the longitudinal modes, namely the so-called stress-diffusion coupling.!' Bi^ The present analysis 
shows that the dynamical asymmetry coupling between the transverse modes is also important for 
the hydrodynamic relaxation process. Previously Edwards and Freed also discussed the same effect 
for polymer solutions by a microscopic arguments .EBBi^ However, they considered a system where 
polymers are assumed not to be entangled. The present result is due to the collective motion on a 
scale of the viscoelastic length, ^„e; resulting from the entanglement of the polymers. 

We now investigate a more general situation. Within the linear regime, the total shear stress 
a\^{t) can be written as 

rt 



J — oo 



where 



G^\t) is called the shear relaxation modulus and fulfllls the following relationl^ 



lUJ 



V0+ r dtG^^\t)e-'^' 
Jo 



With the aid of eqs. (31) and (32), the inverse Laplace transformation of eq. (41) gives 



Gif)(i) = ^exp 
To 



To 



It should be noted that the relaxation time is given not by tq but by tq/[1 + {ivek)"^]- 
We can rewrite eq. (39) in real space as 

rt 



a{r,t) = riok{r,t)+ f dt' f dr'G^P\t - t' ,r - r')k{r' ,t'), 

J —oo J 



(39) 



(40) 



(41) 



(42) 



(43) 



where 



K(r,t) = Vv^(r,t) + [Vv^(r,t)]t, 
and G^P\t — t',r — r') is given by 



(44) 
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To 


exp 


To 





it-t') Ir-rf TO 



To 



(45) 
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Equations (43) and (45) explicitly show that the shear stress o"(r, t) is dominated by the velocity 
fluctuations of the distance |r — r'| ~ ©(Cue) and of the time t — t' ^ 0{tq). If we take the limit 
^ye 0, eq. (43) reduces to the well-known form as follows: 

k{r,t'). (46) 



a{r, t) = riQk{r, t) + — dt' exp 

To J-oo 



it - t') 



To 



In this study, we have investigated the relaxation dynamics of the concentration fluctuations and 
the velocity fluctuations by the linearized two-fluid model equations. The present results can be 
applied to the polymer solutions in a one-phase region far from the critical point. However, in the 
vicinity of the critical point the nonlinear hydrodynamic interaction, resulting from the streaming 
type mode coupling, becomes dominant.^-* Taking such nonlinear hydrodynamic interaction into 
account by means of the mode coupling theory, we will report the critical dynamics of polymer 
solutions elsewhere.® 
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